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THE  DISTRIBUTION  Gf  3UBMUNITION  ARRIVAL  TIMES 


1.  INTRODUCTION 

Recently,  interest  has  been  revived  in  a  comparison  of  the 
antipersonnel  effectiveness  of  an  artillery  shell  equipped  with  a 
mechanical  time  (MT)  fuzed  improved  conventional  munition  (ICM)  versus 
one  which  uses  a  proposed  proximity  (VT)  fuze.  For  such  a  comparison 
one  must  recognize  that  in  reality  an  enemy  target  will  react  to  the 
audible  detection  of  an  incoming  projectile  by  seeking  protective  cover 
in  order  to  reduce  his  exposure  and  therefore  enhance  his  chance  of 
survival.  Thus,  in  addressing  the  effectiveness  of  volley  fire  in  this 
situation  one  must  examine  for  each  candidate  fuze,  subsequent  to 
warning,  both  the  distribution  of  arrival  time  of  submunitions  in  the 
target  area  and  the  target  personnel  reaction  t*^e  distribution  in  order 
to  be  able  to  realistically  assess  the  ability  of  the  target  personnel 
to  achieve  adequate  cover  before  a  significant  number  of  submunitions 
have  functioned. 

This  report  considers  a  theoretical  investigation  of  the  first 
of  these  distributions,  the  arrival  time  of  submunitions  in  the  target 
area.  The  approach  which  is  used  accounts  for  three  primary  sources  of 
variability,  variation  in  time  to  payload  ejection  (due  to  the  fact  that 
volley  rounds  will  not  function  simultaneously),  variation  in  burst 
height  (due  to  the  fact  that  these  rounds  may  function  at  different 
altitudes),  and  finally  variation  in  submunition  time  to  fall  (due  to 
the  fact  that  given  ejection,  all  submunitions  from  a  given  round  will 
also  not  function  simultaneously).  The  assumptions  required  concerning 
the  specific  functional  form  for  the  probability  distributions  are 
based  on  actual  test  data.  The  distribution  of  time  to  payload  ejection 
is  assumed  to  be  a  negative  exponential  which  reasonably  approximates 
the  test  results  of  the  HELBAT  I*  firings  conducted  by  the  Human  Engineer¬ 
ing  Laboratory. 

♦Human  Engineering  Laboratory,  TM  24-70,  "HEL  Battalion  Artillery  Test," 
Sep  1970,  G.  L.  Harley,  D.  J.  Giordano. 
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Similarly,  the  uniform  assumption  for  the  distribution  of  submurition 
time  to  fall  over  the  interval  from  first  to  last  submunition  arrival 
from  a  given  round  does  not  appear  to  be  inconsistent  with  existing 
firing  record  data.*  Details  concerning  the  rationale  used  in  choosing 
specific  parameter  values  can  be  found  in  Section  4.  Section  2  contains 
a  brief  account  of  the  general  methodology  which  includes  a  description 
of  the  basic  underlying  assumptions.  This  is  then  followed,  in  Section 
3,  by  a  detailed  description  of  the  mathematics  involved  in  the  derivation 
along  with  the  final  general  mathematical  expression  for  the  distribution 
in  question.  Finally,  Section  4  is  devoted  to  a  specific  application 
which  compares  the  submunition  arrival  time  distributions  for  several 
munition  configurations  of  current,  interest. 

2 .  GENERAL  ^THODOLOGY 

2.1  Basic  Assumptions.** 

•  Time  is  initialized  to  the  time  of  ejection  of  the  submu¬ 
nitions  from  that  round  of  the  volley  which  is  first  to  function.  The 
noise  signature  associated  with  this  initial  ejection  is  assumed  to  con¬ 
stitute  first  warning.*** 

•  Given  warning,  the  time  to  ejection  for  any  subsequent 
round  is  a  continuous  random  variable  with  probability  density  function 
(p.d.f.)  g3(te)  ,  0  <  te  <  °°.  For  the  results  presented  in  this  report 
this  p.d.f.  is  assumed  to  be  of  the  form 

_  !e 

g3(te)  =  ~  e  °  0  <  te  <  »,  0  >  0. 

•  Burst  height  is  a  continuous  random  variable  with  p.d.f. 

*FR  #  5296,  FR  #  S419  Yuma  Proving  Ground. 

**For  further  amplification  of  the  practical  significance  of  these 
assumptions  the  reader  is  referred  to  Section  4.2. 

***Note  that  no  correction  has  been  made  for  the  time  required  for  the 
sound  of  the  initial  burst  to  reach  the  ground.  If  necessary,  however, 
the  results  of  this  effort  can  easily  be  adjusted  to  reflect  the 
uncertainty  in  when  warning  actually  occurred. 
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gjCh),  a  <  h  <  b.  The  specific  parameter  values  of  this  distribution 
will  depend  on  fuze  type  and  will  also  generally  depend  on  the  angle  of 
fall  of  the  projectile.  Note  that  in  this  report  it  is  assumed  that 
there  exists  a  mean  burst  height  for  the  volley  fire  and  gj(h)  represents 
the  distribution  about  this  mean.  No  attempt  has  been  made  to  address 
the  distribution  of  mean  burst  heights  over  a  series  of  occassions. 

•  Given  submunition  ejection,  submunition  time  to  fall  is  a 
continuous  random  variable  with  conditional  p.d.f.  g->(t-/h),  c(h)  <  t^.  <  d(h). 
For  the  results  presented  in  this  report  this  p.d.f.  is  assumed  to  be 
uniform  where  the  time  to  first  submunition  function,  c(h),  and  the 
spread,  d(h)-c(h),  are  both  linearly  dependent  on  burst  height.  Let 
c(h)  =  0t1  +  Bjh  and  d(h)-c(h)  =  a2+B2h  (“l-- ®»a2>  0,  &2  >_  0).* 

Then  the  p.d.f.  of  submunition  time  to  fall,  is  given  by 


g2(tf/h)  = 


1 

V62h 


Oj+B^h  <  tf/h  <  (61+82)h. 


2 . 2  General  Theory. 

Based  on  the  previous  assumptions  the  time  to  fall  distribution 
can  be  obtained  by  using  the  following  relationship: 


g4(V 


f  g,(tf/h)g. (h)  dh, 
R(tf) 


where  Rft^)  is  the  range  of  integration.  Basically  this  is  the 
expression  for  submunition  time  of  fall  averaged  over  height.  Now  the 
total  time  from  initial  warning  to  submunition  hurst,  T,  is  a  random 
variable  which  can  be  expressed  as  the  sum  of  two  independent  random 
variables,  T  and  That  is,  T  =  T  +T-  where  T  is  the  time  to 
round  function  after  initial  warning  with  p.d.f.  Cte)  and  T^  is  the 
time  to  submunition  burst  with  p.d.f.  g^(t^).  Since  these  variables 

♦Note  that  £or  the  derivations  provided  in  Section  3  only  the  constant 
spread  case  (i.e.,  &2  =  0)  was  considered. 
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are  independent  their  joint  p.d.f.  is  given  by  g6(tg,  tf)  =  83^)  *84^^)  • 
Using  the  change  of  variable  technique  the  joint  p.d.f.  of  T  and  T  , 

87^' *s  ®6^e*  ^~<ce^  by  the  appropriate  Jacobian  of  the 

transformation.  The  marginal  p.d.f.  of  T,  the  total  time,  can  then  be 
obtained  by  integrating  out  the  time  to  round  function  variable.  Thus, 
the  p.d.f.  of  total  submunition  arrival  time  (after  first  ejection)  is 

g5(t)  =  f  g6(t  ,  t-t  )  dt 
5  R(t)  6  e  e  e 

Its  corresponding  cumulative  probability  function, 
t 

F(t)  =  /  g.Cx)  dx, 

0  * 

depicts  the  percent  cf  submunitions  which  have  burst  within  t  seconds 
after  initial  warning. 

The  -following  section  provides  the  details  involved  in 
combining  the  basic  assumptions  of  Section  2.1  with  the  general  theory 
in  order  to  obtain  specific  expressions  fer  those  distributions  of  interest. 


3.  THEORETICAL  DETAILS 

3. 1  The  Distribution  of  Submunition  Time  to  Fall. 

The  probability  density  function  of  submunition  time  to  fall 

is  given  by  the  expression  g. (t .)  =  /  g^Ct^/hJg. (h)  dh.  The  range 

R(tf) 

of  integration,  R(t^.),  is  functionally  dependent  on  t^  and  can  be 
described  explicitly  by  examining  the  composition  of  the  two  dimensional 
space  over  which  the  product,  g2 (t^/h) • g1 (h) ,  takes  on  nen-zero  values. 
This  space  consists  of  all  two  dimensional  points  (t^,  h)  such  that 
both  a  <  h  <  b  and  c^+f^h  <  tf  <  (a1+a2)+  (B^+gph.  Two  distinct  cases 
must  be  considered  in  describing  R(t^).  These  are: 

a.  (a^c^)*  (gj+g^a  <_  cr^+B^b  (Case  A), 

b.  (a1+a2)+(81+62)a  >  o^+Bjb  (Case  B) . 

The  geometry  pertaining  to  these  two  cases  is  pictured  in  Figure  3.1. 
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Figure  3.1  Geometry 


(1) 


3.1.1  Case  A.  For  those  conditions  when  < 

al+3lb’  g4^f^  is  given  by  one  of  the  following  expressions  depending 
on  the  range  of  t^. 


f_°i 


"i 

g2Ctf/h)g1 00  dh  when  Oj^a  <  tf  <  (ai+a2)+  (3i+so)a 

« 


^(tf) 


tf'al 


< 


r 


g2(tf/h)g1  DO  dh  when  (a^c^)* (51+S0)a  <_  tf  <  a  +£  b 


tf~(°l>a2) 

V82 


V 

! 


g2(tf/h}g1(h)  dh  when  Oj+Bjb  <_  tf  <_  (a1+a,)+(£1+/?)b. 


tf-  (a1+a2) 


V.  V62 


Now  for  the  constant  spread  case  (i.e.,  S2  «  0),  g2(tf/h)  =  for  all  points 
(tf,h)  critical  to  the  evaluation  of  the  integrals  defined  in  equation  Cl). 
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3.1.2  Case  B.  Similarly  for  those  conditions  when  (a^*a2)*($  )  a  > 

“l+Sjb,  «4(tf)  is  given  by  one  of  the  following  expressions  depending  on 
the  range  of  t£. 


(0 


I 


g2(tf/h)g1(h)  dh  Qj+Bja  <_  tf  <  a^Bjb 


MV  1  fg2(tf/h)gl(h)  dh  <  t£  <  (o1*o2)+(BI+S2)a 


4 


t 

I 


g2vtf/h)gl(h)  dh  (‘VV+(Wa  <  tf  <  (a1+a2)+(8]+g?)b. 


tf' (“i+V 
81+S2 

For  the  constant  spread  case  (i.e.,  87=0)  equation  (4)  reduces  to 


(5)  MV  H 


1_G  /V°i 
«,  6 


—  when 
a2 


when  ai+eia  i  tf  <  c^+Bjb 


al+8lb  ll£  <  ^ai+V+8la 


when  (a1+a2)+B1a  <_  t£  <  (e^+c^j+B/.;, 


where  G^y)  =  /  g^h)  dh,  a  <_  y  <_  b. 
a 
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3.2  The  Distribution  of  Total  Subaunition  Arrival  Time  (Constant 
Spread  Case:  S  ..=(?). 


3.2. 1 _ Introduction.  As  described  in  Section  2.2  the 

marginal  p.d.f.  of  total  tine  to  fall,  T,  is  given  by: 


(o)  M 


S(t>  -  J  W 


)  dt 
e  e 


R(t) 


/  «3(tc)'84lt-tc)  J£o’ 


R(t) 


where  the  range  of  integration,  R(t),  is  functionally  dependent  on  t  and 

can  be  described  explicitly  for  a  given  value  of  t.  In  order  to  describe 

R(t)  it  is  helpful  to  examine  the  two-dimensional  region  over  which  the 

joint  probability  of  (T,  Te),  g_ 'tQ) *g4 (t-tg) ,  takes  on  non-zero  values. 

Noting  that,  for  the  constant  spread  case  (i.c.,  t.=0),  g_(t  )  >  0  for 

0  1  ce  <  °°>  fij  (*f)  >  0  for  +  <_  tf  <_  jb,  and  T  =  ff+Te» 

it  follows  that  the  region  in  question  consists  of  all  points,  (t,  tg), 

such  that  simultaneously  0  tp  <  «  and  ja  +  tg  t  <_  (aj+a7)+Bjb+te. 

This  region  is  depicted  in  Figure  3.2  and  consists  of  all  points  in 

the  shaded  region.  Using  this  figure  it  can  be  seen  that  the  range  of 

integration  for  t=t  is  given  by:  R(t  )={t:0<t  <t 

Similarly,  for  t  =  t^,  the  range  of  integration  is  R(tj)  =  {tp:  ti-aj-a2“ 

B^b  <_  t  < t^-a^-B^a} .  Before  the  integration  involved  in  equation  (6) 

can  be  performed  the  region  specified  in  Figure  3.2  must  be  further 

partitioned  into  subregions  corresponding  to  the  various  definitions  of 

g4(t^.J.  As  in  Section  3.1  both  Case  A  and  Case  B  must  be  considered. 

These  subregions  are  shown  in  Figure  3.3.  For  Case  A  -  Region  I, 

a.+B.a+t  <  t  <  a.+u.+B.a+t  and  thus  a.+B.a  <  t-t  <•  cx  +a.+6.a. 
lie—  -121  e  11—  e  121 

Now  using  the  fact  that  t^.  =  t-te  and  the  definition  of  g^Ct^.)  contained 
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Figure  3.2  Points  of  Positive  Probability. 


2  /*'  uj  \ 

in  equation  (5)  it  follows  that  g,(t-t  )  =  — G,| — - - )  in  Region 

4  e  a2  1\  J 

I  -  Case  A.  A  similar  argument  car:  be  used  for  each  of  the  five 
remaining  special  cases.  The  results  of  each  of  these  investigations 
is  contained  in  Table  3.2. 

Having  obtained  the  above  partition  one  can  iow  proceed  with 
the  integration  involved  in  equation  (6).  Cases  A  an  1  B  will  be 
considered  separately. 

3.2.2  Case  A.  The  integration  in  question  is  that  involved 
in  the  following  expression: 


g5«) 


f 


R(t) 


g.  (t  ) • g, ft-t  )  dt 
e  ■  e  ( 


For  the  sake  of  clarity  this  integration  will  be  performed  in  stages, 
each  stage  representing  the  integration  over  a  specific  range  of  t. 

In  all  cases  Figure  3.4  anc  Table  3.2  will  be  useful  gui.'es. 

For  a^Sja  1  t  <  V  W* 


g5(t) 


t-  ( ij+^a) 

)  83(t0)84(t-te)  <lte. 

0 


Now  for  the  specific  distributional  assumption  made  in  Section  2.1  this 
expression  reduces  to: 


t-  Osj+Bja) 


0 
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TABLE 


1 

] 
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Similarly,  for  a  +a,+  p.  a 


This  expression  reduces  to 


Thus,  for  this  range  of  t. 


23 


Note  that  each  of  the  integrals  involved  in  the  above  representation 


of  g5(t)  is  of  the  form 


(7) 


s 


a2 


G, (m*nt  )  dt  . 
1  e'  e 


r 


Consider,  for  example,  the  first  integral  involved  in  the  expression  which  is 


0 


1  1 

In  this  case,  r  =  0,  s  =  t-(a  +3.a),  m  =  —= — ,  and  n  =  -  — -  . 

it  dj  ^ 

Using  the  technique  of  integration  by  parts  the  following  general  solution 
for  (7)  results: 


s 


m+nr 


For  the  special  case,  (8),  it  follows  that 


t-Caj+gja) 


C1 


Similarly,  each  of  the  integrals  involved  in  the  previous  expression  for 
g<_(t)  can  be  reduced  and  after  some  algebraic  simplification  they  form 
the  following  expression  for  g5(t). 

For  a1+“ia  <  t  <  (Oj+a^+Sja, 


a 
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Similarly , 


5.2.3  Case  B.  By  repeating  the  methodology  presented  in  Section  5.2.2, 
using  the  geometry  and  functional  definitions  appropriate  for  Case  B,  the 
following  formulation  for  g-(t)  results: 


27 


\ 


a  a 


4.  SPECIFIC  APPLICATION 
4. 1  Introduction. 

The  following  is  a  description  of  a  specific  application  of  the 
theory  described  in  Sections  2  and  3.  Note  that  where  possible  existing 
data  were  used  to  construct  the  underlying  distributions.  As  further  data 
become  available  vhe  results  can  easily  be  updated  to  reflect  this  additional 
information.  In  any  case,  this  description  should  provide  the  reader  with 
a  comparison  of  the  distribution  of  submunition  arrival  times  for  several 
specific  munition  configurations  of  current  interest,  a  clearer  understanding 
of  how  the  methodology  might  be  applied,  and  some  indication  of  the 
sensitivity  of  the  results  to  .'ariations  in  the  basic  parameter 
values. 

The  application  considers  a  battery  volley  of  artillery  fire  using 
Improved  Conventional  Munitions.  A  comparison  is  made  of  the  per  cent 
of  submunition  bursts  as  a  function  of  time  after  initial  warning  for 
several  specific  munition  configurations.  As  specified  in  Section  2 
initial  warning  occurs  (i.e.,  time  is  initialized)  at  the  time  of 
ejection  of  the  submunitions  from  that  round  of  the  volley  which  is 
first  to  function.*  Three  munition  configurations  were  considered  for 
the  ISSmm,  M449A1  ICM  with  specific  characteristics  as  listed  in  Table  4.1. 

‘Note  that  no  correction  has  been  made  for  the  time  required  for  the  sound 
of  the  initial  burst  to  reach  the  ground.  If  the  uncertainty  associated 
with  actual  warning  is  thought  to  be  significant  the  warning  time  axis 
should  be  adjusted  appropriately. 
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TABLE  4.1  MUNITION  CONFIGURATIONS 


Configuration 

Fuze 

Angle  of  Fall 

I 

Mechanical 

Time  (MT) 

XM  577 

30° 

II 

MT 

XM  577 

50° 

III 

Proximity  (VT) 

Not 

Applicable 

The  distributional  assumptions  required  by  the  theory  of  Sections  2  and 
3  will  now  be  discussed  for  each  of  these  cases. 


4.2  Distributional  Assumptions. 

4.2.1  g,(tc)  ~  Round  Arrival  Time.  In  practice  there  are 
many  sources  of  error  which  contribute  to  the  variation  in  round  arrival 
time*  when  firing  a  battery  volley.  The  major  sources  of  variability 
are  the  inability  of  the  firing  battery  to  conduct  a  synchronized  firing 
and  delivery  error  sources  such  as  tube-to-tube  velocity  variations, 
velocity  variations  within  a  given  tube,  aiming  variations,  etc. 

The  contribution  of  differences  in  fuze  type,  on  the  other  h.'nd,  will 
most  likely  be  relatively  insignificant  compared  to  those  of  the  factors 
mentioned  above.  Thus,  for  the  purposes  of  this  example  Ctg)  is 
assumed  to  be  independent  of  fuze  type.** 

As  stated  in  Section  2.1  the  general  form  assumed  for  the 
p.d.f.  of  round  arrival  given  warning  is 


0  <  t  <  0 

e 


♦Round  arrival  time  is  considered  to  be  the  time  to  round  function 
subsequent  to  the  function  of  the  first  round  of  the  volley. 

**Note  that  this  assumption  was  made  by  choice  and  should  not  be  construed 
as  a  restriction  inherent  in  the  methodology. 
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For  each  of  the  cases  specified  in  Table  4.1,  results  are 
obtained  for  two  values  of  the  parameter  9  (0  =  1  and  9=2).  In  the 
case  where  8=2,  gjftg)  reasonably  approximates  the  round  arrival  time 
data  obtained  from  HElBAT  I.*  On  the  other  hand,  certain  situations  may 
exist  where  a  slightly  increased  likelihood  of  shorter  arrival  times  is 
more  realistic.  6=1  reflects  such  a  situation.  In  any  case, 
consideration  of  both  values  provides  insight  into  the  sensitivity  of 
the  methodology  to  this  parameter.  The  practical  significance  of  the 
values  6=1  and  0=2  can  be  better  appreciated  by  examining  the 
cumulative  distribution  functions**  displayed  in  Figure  4.1.  Examining 
the  curves  it  is  seen  that  90%  of  the  rounds  from  the  voliey  will  function 
within  a  span  of  2.3  seconds  when  6=1  and  within  a  span  of  4.65  seconds 
when  5=2.  As  stated  earlier,  if  other  values  of  0  appear  to  be  more 
realistic  as  further  information  is  obtained  then  these  cases  should 
also  be  examined. 

4.2.2  g^Ch)  -  Burst  Height.  Projectile  burst  height  is 
assumed  to  be  a  normal  random  variable  with  p.d.f. 


(h-p)2 

2o2 

g  (h)  =  — r-  e 
1  /27  a 


2 

-«■  <  h  <  00 ;  _re  <  p  •*  to;  c“  >  0. 


*Human  Engineering  Laboratory,  TM  24-70,  "Human  Engineering  Laboratories 
Battalion  Artillery  Test,"  Sep  1970,  Gary  L.  Harley,  Dominick  J. 
Giordano. 


~F(to) 


f  -te/e 

J  g3('te')  (lte  =  1-e 


6>0;0<t  <°° 

o 
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In  practice  burst  heights  are  limited  in  that  rhey  will  neither  occur 

below  the  surface  (i.e.,  h  >_  0)  nor  above  some  limit,  say  L,  determined 

by  the  ballistic  altitude  capabilities  of  the  gun/shell  system  (i.e., 

h  <_  L).  Thus  the  normal  p.d.f.,  with  infinite  range,  is  not  a  precise 

representation  of  the  true  burst  height  distribution.  In  all  cases 

considered,  however,  the  normal  approximation  is  adequate  since  points 

a  » 

a  and  b  (a  <  b;  a  >  0,  b  <  L)  do  exist  such  that  I  gj(h)  dh  and  /  gj  (h)  dh 
are  both  approximately  equal  to  zero. 


~> 

The  form  of  the  distribution  (i.e.,  the  values  of  u  and  o  ) 
will  generally  depend  on  the  various  components  of  the  delivery  error 
(precision,  MPI,  fuze),  the  fuze  type,  and  the  angle  of  fall  of  the 
projectile.  Specific  values  for  the  three  munition  configurations 
considered  here  were  taken  from  a  recent  report  by  the  USAMJCOM  ORG* 
and  are  listed  in  Table  4.2.  Details  concerning  the  rationale  supporting 
these  values  can  be  found  in  the  referenced  document. 


TABLE  4.2  NORMAL  PARAMETERS 


Configuration 

Mean  Burst  Height 
(u) 

Standard  Deviation 
(o) 

I 

23S  meters 

25.7  meters 

II 

350  meters 

76.3  meters 

III 

187  meters 

4.8  meters 

4d-.L.J.2. 

(c-/h)  -  Submunition  Time 

to  Fall.  As  specified  in 

Section  2.1,  submunition  time  to  fall  (given  round  function)  is  assumed 
to  be  a  uniform  random  variable  where  the  time  to  first  submunition 


function,  c(h),  and  the  spread,  d(h)-c(h),  are  both  linearly  dependent 
on  the  height  that  the  round  functions.  Specifically,  the  methodology 
of  Section  2.3  uses  c(h)  =  c^+B^h  and  d(h)-c(h)  =  (constant  -  independ¬ 
ent  of  burst  height).  For  this  example  both  a  six  second  (a^=6)  and 
a  four  second  (c^^)  spread  were  considered.  In  addition  results  were 

*"A  Comparison  of  the  Effectiveness  of  Improved  Conventional  Munitions 
Equipped  with  Mechanical  Time  and  Proximity  Fuzes  (U)",  DRAFT,  USAMUCOM 
ORG,  November  1972,  CONFIDENTIAL. 
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obtained  for  two  different  relationships  for  time  to  first  submunition 
burst.  These  are  c(h)  =  .P23h  and  c(h)  =  .018h  where  the  time  to  first 
burst,  c(h),  is  expressed  in  seconds  and  burst  height,  h,  in  meters. 

The  .023  relationship  was  chosen  since  it  reasonably  approximates 
existing  firing  record  data*  for  the  M43  type  bomblet.  The  .018 
relationship,  on  the  other  hand,  was  examined  since  it  reflects  the 
rate  of  fall  essumed  in  the  previously  referenced  USAMUCOM  ORG  study  and 
seemed  to  be  a  reasonable  alternative  to  consider  in  examining  the 
sensitivity  of  the  methodology  to  this  assumption. 

4.3  Results. 

The  methodology  of  Section  3.2  has  been  exercised  for  eight 
specific  cases.  Each  case  represents  a  comparison  of  the  distribution 
of  submunition  time  to  fall  for  three  configurations  of  the  155mm,  M449  A1 
ICM.  Configurations  i  and  II  involve  the  XM577  mechanical  time  fuze  at  a 
30°  and  S0°  angle  of  fall,  respectively,  while  configuration  III  involves 
a  proposed  proximity  fuze.  The  distributional  parameter  values  considered 
in  each  of  the  eight  cases  are  listed  in  Table  4.3. 

TABLE  4.3  SPECIFIC  CASES 


Case 

Time  to  Ejection 
Parameter  (6) 

Spread  From  First  to 

Last  Submunition  Function  (a^) 

Height  Dependence 
*1 

1 

2 

6  seconds 

.023 

9 

2 

4  seconds 

.023 

3 

1 

6  seconds 

.023 

4 

1 

4  seconds 

.023 

5 

2 

6  seconds 

.018 

6 

2 

4  seconds 

.018 

7 

1 

6  seconds 

.018 

8 

1 

4  seconds 

.018 

The  rationale  for  the  selection  of  these  values  of  9,  cs,,  and  is  contained 


*FR  #5296,  FR  #5419,  Yuma  Proving  Ground. 
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in  the  discussion  of  Section  4.2.  Note  that  values  other  than  the  nomi- 
nal  values  have  also  been  examined  in  order  to  get  soae  indication  of 
the  sensitivity  of  the  results  to  variation  in  the  basic  parameter  values. 
Figures  4.2  thru  4.9  display  the  results  of  this  evaluation.  In  each 
case  the  cumulative  probability  distribution*  corresponding  to  the  p.d.f. 

of  total  time  to  fall,  g5(t),  is  provided.  It  is  seen,  for  Case  1  for 
example,  from  Figure  4.2  that  for  the  initial  conditions  specified  in 
Table  4.3  only  35  per  cent  (Munition  Configuration  III)  and  20.5  per 
cent  (Munition  Configuration  I)  of  the  volley  submunitions  will  burst 
within  the  first  eight  seconds  after  initial  warning  (first  round 
submunition  ejection).  In  fact  in  all  cases  there  wii.'.  be  considerable 
time  (4-7  seconds)  for  personnel  to  react  before  many  of  the  submunitions 
from  the  volley  detonate.  Additionally  this  estimate  may  even  be 
conservative  since  in  the  operaticnal  environment  there  mav  possibly 
be  an  additional  2  to  3  seconds  of  warning  due  to  the  incoming  whistling 
of  the  round  prior  to  submunition  ejection. 

Another  observation  that  can  be  made  is  that  the  results  are 
quite  sensitive  to  the  parameters  0,  and  Bj.  Thus  the  variability 
in  round  arrival  time  and  variation  in  time  to  submunition  impact  do 
significantly  affect  the  ultimate  distribution  of  submunition  arrival 
time.  Failure  to  recognize  and  account  for  these  sources  of  variability 
could  seriously  bias  the  results  of  any  study  evaluating  weapons  effects. 

5.  SUMMARY 

In  assessing  the  effectiveness  of  an  artillery  battery  volley 
using  Improved  Conventional  Munitions  in  a  situation  where  the  personnel 
in  the  target  area  may  react  to  seek  protective  cover,  one  must  consider 
both  the  distribution  of  arrival  time  of  submunitions  in  the  target  area 
and  the  reaction  time  distribution  for  the  target  personnel.  The  first 

-  t 

o 

*The  cumulative  distribution  function  F(t  )  =  f  g_(t)  dt  depicts  the  per 

0  o  0  ... 

cent  of  submunitions  which  have  burst  within  t  seconds  after  initial 
warning. 
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of  these  has  been  addressed  in  this  report  and  a  mathematical  expression 
for  the  distribution  of  submunition  arrival  time  is  provided  in  Section  3. 

In  the  subsequent  analysis,  comparing  several  munition  configurations  of 
current  interest,  in  Section  4  it  is  shown  that  variability  in  time  to 
round  function  and  variation  in  time  to  submunition  fall  do  significantly 
affect  the  distribution  of  submunition  burst  time.  Further  even  when 
employing  the  VT  fuze  there  will  in  most  cases  be  considerable  time  for 
personnel  to  react  before  many  of  the  submunitions  from  the  volley  detonate. 

To  realistically  assess  the  degree  of  surprise  which  can  be 
achieved  with  alternative  munition  configurations  results  such  as  those 
provided  should  be  coupled  with  information  pertaining  to  the  ability 
of  target  personnel  to  react  to  audible  warning  in  seeking  protective 
cover.  However,  the  magnitude  of  the  warning  times  alone  make  one 
seriously  question  whether  any  significant  degree  of  surprise  can  be 
achieved  by  using  any  of  the  munition  configurations  considered  in  this 
report.  This  opinion  is  reinforced  by  the  small  quantity  of  test  data 
available  in  HEL  Technical  Notes  7-68  and  1-71  where  troop  posture 
sequences  were  tested  for  dismounted  'armored)  infantrymen  and  for 
artilleiy  crews  respectively.  These  reports  both  indicated  that  with 
4-6  seconds  of  time  to  react  target  personnel  would  be  able  to 
significantly  decrease  their  total  exposed  presented  area.  As  further 
data  become  available  on  target  personnel  posture  distributions  an 
additional  analysis  may  be  required  to  completely  summarize  the  ICM 
effectiveness  envelope  for  area  personnel  targets. 
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Figure  4.2  Cumulative  Distribution  of  Submunition  Arrival  Times. 


Figure  4.3  Cumulative  Distribution  of  Submunition  Arrival  Times. 


Figure  4.4  Cumulative  Distribution  of  Submunition  Arrival  Times. 


Figure  4.5  Cumulative  Distribution  of  Submunition  Arrival  Times, 


Figure  4.6  Cumulative  Distribution  of  Submunition  Arrival  Times. 


Figure  4.7  Cumulative  Distribution  of  Submunition  Arrival  Times. 


Cumulative  Distribution  of  Submunition  Arrival  Times. 


Figure  4.9  Cumulative  Distribution  of  Submunition  Arrival  Times. 


